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CN : Abstract 

U 

^^ , The reaction-diffusion system of Fitzhugli Nagumo is considered. The initial- 

^ ' boundary problems with Neumann and Dirichlet conditions are analyzed. By 

^>0 , means of an equivalent scmilincar intcgrodiffcrential equation which character- 

izes several dissipative models of viscoelasticity, biology, and superconductivity, 
some results on existence, uniqueness and a priori estimates are deduced both 
r^ ' in the linear case and in the non linear one. 
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^ ! 1 Introduction 

'NT 

^, 

'nT . As it is well known, the FitzHugh - Nagumo model (FHN) has been introduced 

ir I by FitzHugh (1961) and Nagumo et al.(1962) as a simplified description of the 

excitation and propagation of nerve impulses [1,2]. 

Of course, there are many other important physiological applications in addition 
to that for the propagation of nerve action potentials modeled by the FHN system. 
One such important application is related to the waves that arise in muscle tissue, 
^_' particularly heart muscle and particularly interesting is the biophysical phenomenon 

of reentry which occurs in the excitable cells of the heart. [3] Another example is the 
reverberating cortical depression waves in the brain cortex. [4] As for the analysis of 
the (FHN) model, there exists a large bibliography both in the linear case and non 
linear problem. [5,7-9]. 

The model is given by the set of p.d.e [4, 10]: 
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(1.1) 

I rli) 

(3v. 

where u{x,t) models the transmembrane voltage of a nerve axon at distance x and 
time t, V {x,t) is an auxiliary variable that acts as recovery variable. Further, the 
diffusion coefficient e and the parameters b, (3 are all non negative [11,12]. 
Besides, the function f{u) has the qualitative form of a cubic polynomial 

(1.2) f{u) = — au + (p{u) with ip = u'^ {a + 1 — u) 

where the parameter a is the threshold constant and generally one has < a < 1. 
However, in many papers phenomena for the (FHN) system has been investigated 
considering a as a non homogeneous function and a < 0. [13, 14] 

Moreover, Denoting by vq the initial value of v, system (1.1) can be given the 

form of the following integro differential equation 

(1.3) Cu = ut — euxx + au + b / e~^^^~'^' u{x, r) dr = F{ x, t, u (re, t) ) 

Jo 

as soon as one puts: 

(1.4) F{x,t,u) = ip{u) - vo{x)e~'^K 

Equation (1.3) describes the evolution of several physical models as motions of 
viscoelastic fluids or solids [15-17]; heat conduction at low temperature [18], sound 
propagation in viscous gases [19]. 

Moreover, it occurs also in superconductivity to describe the Josephson tunnel 
effects in junctions. In this case the unknown u denotes the difference between the 
phases of the wave functions of the two superconductors and it results: 

(1.5) euxxt - utt + Uxx - OiUt = sinu +7 

where 7 is a constant forcing term that is proportional to a bias current, and the 
e — term and the a — term account for the dissipative normal electron current flow 
along and across the junction, respectively [20]. 



From (1.3) one obtains the equation (1.5) as soon as one assumes 

1 a ^ 1 

1.6 a = a -- b= -- (3 = - 

e e e 

and F is such that 

(1.7) F{x,t,u) = - / e"^(*"^) [7 + senn(x,r) ] dr. 

Jo 

As (1.6) show, in the superconductive case the constants a, h could be negative 
too. 

In this paper Neumann and the Dirichlet problem for (1.3) are considered. By 
means of the fundamental solution KO(x,t) of the operator L, which has already 
been determined in [21], in the linear case, the explicit solutions are obtained. When 
the source term F is non linear, an appropriate analysis of the integro differential 
equation implies results on the existence and uniqueness of solutions. Besides, in 
both the linear case and the non linear problem, a priori estimates of the solutions 
are achieved. The results are applied to (FHN) system. 

2 Statement of the problems and Laplace transform 

If T is an arbitrary positive constant and 

Or={(x,t):0<x<L; 0<t<r}, 

let (-Pat) the following Neumann initial - boundary value problem related to equation 
(1.3): 

ut — euxx + au + b Jq e~^^^~'^' u{x, r) dr = F{x, t, u) {x, t) E VLt 
(2.1) { 'u(x,0) =uq{x) X G [0,L], 

'Ux(0,t) = ^i{t) Ux{L,t) = Mt) 0<t<T. 



In excitable systems this problem occurs when two-species reaction diffusion system 
is subjected to flux boundary condition [4]. The same conditions are present in case 
of pacemakers [2] . Neumann conditions are applied also to study distributed FHN 
system [22]. 

Besides, if the Dirichlet conditions are applied, the following initial boundary value 
problem (Pd) holds: 



(2.2) 



Cu = F{x,t,u) {x,t) G fir 

u{x,0) = uo{x) X G [0, L], 

n(0,i) = gi{t) u{L,t) = g2{t) 0<t<T. 



For example, in mathematical biology, those boundary conditions occur when the 
behavior of a single dendrite has to determine and the voltage level has to be fixed. [2] 

When F = f{x,t) is a linear function, problems (Pn), {Pd) can be solved by 
Laplace transform with respect to t. 
If 



(2.3) 



u{x, s) = f^ e ^* u 



g ,. ^{x,t)dt f{x,s) = / ~ e "^ f{x,t)dt, 
^^{s)=!^ e-''^i{t)dt (i = l,2), 



one deduces the following transform {Pn) problem: 



u. 



(2.4) 



l[f{x,s) + uo{x)] 



UxiO,s) = 1^1 (s) Ux{L,s) = ■4^2is), 



where a = s + a -\- -. Letting a = a /e, and considering the followine 

S + /3 to / , 

function 



(2.5) 



%{y,^) 



cosh [a {L — y)] 
2 e G sinh ( cr L ) 
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1^ o 



e ^ 
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n=l 



2nL+y 

e ^ 



2nL — y 

+ e 7^ 



the formal solutions u{x, s) of the problems (Pn) can be given the form: 



u{x,s)= [eo{\x-^\,s)+eo{\x + C\,s)][uo{a + f{C,s)]dC- 
(2.6) Jo 

- 2 e V'l (s) Oo{x, s) + 2e^2 {s) Oq {L-x,s ). 



Analogously, for problem {Pd), one obtains: 



u{x,s)= [ [0o{x + ^,s)- eo{\x-^\,s)][uo{a + f{^,s)]d^ 
(2.7) -^0 

- 2 e gi{s) ^- eo{x, s) + 2 e iIj2{s) ^— 9q{L - x,s). 
ox ox 



3 Fundamental solution Kq{x^ t) and theta function 9q{x^ t) 

The fundamental solution Kq{x, t) of the linear operator C defined in (1.3) has been 
already obtained in [21] and it results: 



J2 



(3.1) Ko{r,t) 



2^ 



Vb ,, e-^'^'~yh,{2^by{t-y))dy 



Vi Jo Vt-y 



where r = |x| /^/e and Jn{z) is the Bessel function of first kind with 

/■oo —r<7 

(3.2) Ct Kq = / e~'* Kq (r, t) dt = — -^- ^es> max{ -a, -13) 
Jo 2^ea 

Indicating by a; = min{a, /3) and 



p-l3t _ -at -^ ^ n-\- B 

(3.3) E(t)= ^_^ >0. ,„=- + .VS^. 

it results [21]: 



rt 



(3.4) \Ko\ < -^7^ [e-^* + btE{t)]; [ dr [ \Ko{x - ^,t)\ d^ < /3o, 



(3.5) / \Ko{x-C,t)\d^ < e-"^ + Vbirt 



e-^K 



In order to obtain the inverse formulae for (2.7), let apply (3.2) to (2.5). Then 
one deduces the following function similar to theta functions: 



(3.6) Oo{x,t) = Ko{x,t) + Y^ [Ko{x +2nL, t) + Ko{x-2nL, t)] 

n=l 

oo 

= ^ Ko{x +2nL, t). 



As consequence, by (2.7), the explicit solution of the /inear problem (Pn) where 
F = f{x,t) is: 



(3.7) u{x,t)= [ [eo{\x-Clt)+eo{x + C,t)]uo{OdC + 

Jo 



2e / Bo{x, t-r) ipiir) dr + 2 e 9o{L-x,t-T) ^2{t) dr 
Jo Jo 

+ I dr f [eQ{\x-ilt-T)+eo{x + i,t-T)] fU,T) dt 

JO JO 



In an analogous way, a similar formula for the problem {Pd) can be obtained: 



(3.8) u{x,t)= [ [9oix + C,t) -9o{\x-C\,t)]uoiC)dC + 

Jo 



-2e / —eoix,t- t) (71 (r) dr + 2 e —eo{L-x,t-T) g^ir) dr 
jQ ox Jo ox 

+ [ dT [ [Ooix + ^, t - t) - eo{\x - ^\, t - t)] f{^,T) dC 
Jo Jo 

Owing to the basic properties of KQ{x,t), it is easy to deduce the following 
theorems: 

Teorema 3.1. When the linear source f{x,t) is continuous in JIt and the initial 
boundary data uo{x), ipi{t) {i = 1,2) [ g-i {i = 1,2)/ are continuous, then problem 
{Pn) [{Pd)] admits a unique regular solution u{x,t) given by (3.7)[(3.8)]. I 



As consequence of the properties of fundamental solution KQ{x,t), various esti- 
mates for w, ut, Ux-.. could be obtained. 

As an example, let evaluate the asymptotic properties of the terms caused by 
the initial datum uo{x) and the source f{x,t). If 

||'Uo||= sup 1^0(3;) I, 11/11 = sup I /(x, t) I, 

0<x<L Qt 

it results: 

Teorema 3.2. When ipi = {i = 1,2), [gi = {i = 1,2)] the solution (3.7)[(3.8)] 
of (P/v),/ (Pd)] for large t, verifies the following estimate: 



(3.9) |n(3;,i)| < 2 [ 11/11 /3o + \\uo\\ {1 + Vbirt) e" 



'Ult 



where a; = min (a, j3) and Po is given by Po = - + vrvft — — —j- 

2(ap)'^/^ 

Proof: Properties of Ko(x,t) imply that: 



(3.10) 



00 ^i 



/ eQ{\x-i\,t)di < V / \K^{\x-i + 2nLlt)\di = 

•^0 n=-oo -^0 

°° /•x+2nL r 

= V / \Ko{y,t)\dy < / \Ko{y,t)\dy. 

„=_^ A+{2n-l)L JSR 

So, applying properties (3.4)2 and (3.5) to (3.7) [(3.8)], the estimate (3.9) follows. I 

4 The FitzHugh-Nagumo model. A priori estimates 

Linear case - If the reaction kinetics of the model can be outlined by means of 
piecewise linear approximations, then one has: 

(4.1) f{u)= 7j(u-a) -u (0 < a < 1) 

where r] denotes the unit- step function [23,24]. Like in (1.2), the linear approxima- 
tion (4.1) involve a linear term — u. As consequence, denoting by fj the constant 
that holds zero or one, the FHN model can be given the form: 



(4.2) 



Ut - eUxx+ U+ V = 7] 

{x,t) G 3? 
^vt+f3v-bu = , 



and estimate (3.9) can be applied. 

Non linear case - Consider now the non linear case defined by (1.2). By means of 
the previous results we are able to obtain integral equations for the two components 
(m, v) in terms of the data. All this implies, also in this case, a qualitative analysis 
of the solution together with a priori estimates. 

At first let us observe that by (1.1)2 one has: 

(4.3) V = foe"^* + b [ e-^(*-^)n(x,r)dr 

Jo 

and this formula, together with (1.4) require the presence of the following convolu- 
tions: 

(4.4) Ki{r,t) = [ e-^^'-^^ Ki.i{x,T)dT (i = l,2) 

Jo 

which explicitly are given by [21]: 

ft e-5^-'^?/-/3(*-2/) / [t^^y-^ , 

(4.5) K,= / -^ U/^ J._i(2V6y(t-y))dy (i = l,2). 

As consequence, together with ^o defined by (3.6), the other two Q functions 

00 

(4.6) 6'i(x,t) = ^ Ki{x^2nL,t) (i = l,2) 

n=— 00 

must be considered. 

To allow a plainer reading let's set 



(4.7) Gi{x,^,t) = 9a\x-^\,t) + ei{x + ^,t) (i = 0,l,2) 
In this manner, owing to (3.7) one has: 

(4.8) u{x,t)= [ [Go(,x,tt) uoiO - Gi{x,C,t)vo{On + 

Jo 

-2e / eo{x,t-T) 7/;i(t) dr + 2 e i 9o{L-x,t-T) ip2{r) dr 
Jo Jo 

+ [ dr [ Go{x,(,t-T) ip[C,T,uitT)] dC. 
Jo Jo 

As for the v component, by (4.3) one deduces: 

(4.9) v{x,t) = voe-^' +b [ [Gi{x,^,t) uoiO -G2{x,^,t)vo{OH + 

Jo 

-2b e / ei{x,t-T) V'i(t) dr + 2be 9i{L-x,t-T) V2(t) dr 
Jo Jo 

+ b [ dr [ Gi(x,e, t-r) (^[C,T,n(e,T)] de- 
Jo Jo 



Let us observe that the kernels Ki(x,t) and K2{x,t) have the same properties 
of ^0(2;,*). In fact [21]: 

Teorema 4.1. For all the positive constants a, be, (3 it results: 

(4.10) / \Ki\ di < E{t)] [ dr [ \Ki\ d^ < (3i 
J^ Jo Jn 

(4.11) / \K2{x-^,t)\ di< I e-^y-^^'-y^ {t-y)dy < tE{t) 
Jn Jo 

where E{t) is defined in (3.3)i and Pi = (a/3)~^. I 



Now, let II z II = sup \z {x, t)\, and let Bt denote the Banach space 
(4.12) Bt = {z{x,t) : z £C{nT), H^H < oo }. 



By means of standard methods related to integral equations and owing to basic 
properties of Ki, Gi (i = 0,1,2) and 9j(n), it is easy to prove that the mapping 
defined by (4.8) [(4.9)] ia a contraction of Bt in Bt and so it admits an unique fixed 
point u{x,t) £ Bt [25,26]. Hence 

Teorema 4.2. When the initial data {uo,vq) are continuous functions, then the 
Neumann [Dirichlet] problem related to the non linear (FHN) system (1.1), (1.2) 
has a unique solution in the space of solutions which are regular in Qt- ' 



Continuous dependence for the solution of (Pn) [(-Pd)] is an obvious consequence 
of the previous estimates. 

As for asymptotic properties, it is well known that the (FHN) system admits 
arbitrary large invariant rectangles S containing (0,0) so that the solution {u,v), 
for all times t > 0, lies in the interior of S when the initial data {uo,Vo) belong to 
S. [?] For this, considering, for example, the case V'l = ^^2 = and letting 

II (/9 II = sup \^p{x, t, u)\, 

then by means of (4.8), (4.9) and owing to the estimates (4.12), (??), (4.10), (4.11), 
the following theorem can be stated: 

Teorema 4.3. For regular solution (n, v) of the (FHN) model, when ^i = "02 = 0, 
the following estimates hold: 



\u\ < 2 [ ||no|| (1 + 7r\/6t) e""^* + ||i;o|| E{t) + /3o \\f\ 
(4.13) 

\v\ < \\vo\\ e-P^ + 2[6(||uo|| + t \\vQ\\)E{t) + 6/3i \\^\\] 



Therefore, when t is large, the effect due to the initial disturbances {uq,vq) is 
exponentially vanishing while the effect of the non linear source is bounded for all t. 

Asymptotic effects of boundary perturbations can be found in [6]. 

10 



All the previous results can be applied to the boundary Dirichlet, too. 
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